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ABSTRACT 


The GISS general circulation model is used to compute global 
monthly mean forecasts for January 1973, 1974, and 1975 from initial 
conditions on the first day of each month, with ocean surface fluxes 
based on climatological mean January sea-surface temperatures. Fore- 
casts are evaluated in terms of global and hemispheric energetics, 
zonally-averaged meridional and vertical profiles, for^^st error 
statistics, and monthly mean synoptic fields. Although it generates 
a realistic mean meridional structure for the month of January, the 
model does not adequately reproduce the observed interannual variations 
in the large-scale monthly mean energetics and zonally-averaged circula- 
tion. The model exhibits no general skill in predicting the monthly— 
mean sea-level pressure field, but it does simulate observed -changes 
in the intensity of the Icelandic low from year to year. For each 
January the model produces a prognostic monthly mean 500-mb height 
field that is superior to climatology and persistence. 

The impact of temporal sea-surface temperature variations on 
monthly mean global forecasts with the GISS model is investigated by 
comparing two parallel forecasts for January 1974, one using climato- 
logical ocean temperatures for the surface flux computations, and the 
other observed daily ocean temperatures. In the one case studied, the 
use of daily-updated sea-surface temperatures produced no discernable 
.beneficial effect on the forecasts, and the total impact on the large- 
scale pressure and temperature fields was small. 
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MONTHLY MEAN FORECAST EXPERIMEOTS WITH THE GISS MODEL 
1. INTRODUCTION 

All dynamical weather predictions presently appear to exhibit a 

similar decay of skill when verified against daily synoptic fields, 

generally showing little or no superiority over climatology after about 
4 

3 to 7 days. (See, e.g,, Miyakoda et al,, 1969, Baumhefner, 1970; 
Miyakoda et al,, 1972; Druyan, 1974; Spar and Atlas, 1975; Druyan et al., 
1975), It has also been shown in a number of theoretical and experi- 
mental studies that inherent uncertainties in the initial state of- the 
atmosphere apparently limit the maximum range of deterministic pre- 
dictability to 2 to 3 weeks, (See, e.g., Smagorinsky, 1969; Lorenz, 1973; 
Leith, 1974), Nevertheless, despite the limited skill of daily prognostic 
maps beyond the short range, it is possible that time-averaged numerical 
predictions may retain some useful skill over longer periods of time. 

There is some- empirical evidence (e.g,, Namias, 1953; 1964) that cer- - 
tain large-scale atmospheric anomalies may preserve their identity over 
periods of a month or even a season, and that time -aver aging, by filter- 
ing out the smaller-scale and shorter-lived components of the atmosphere 
and by diminishing the effects of phase errors, maj' thus yield some de- 
tecjiable skill over extended and long ranges. 

The availability of global general circulation models (GCM's) which 
can be integrated over effectively unlimited time and are theoretically 
unconstrained by lateral boundary errors provides a tempting opportunity 
to test .this hypothesis. In this study we have computed .a set of ex- 
perimental global monthly mean forecasts for three winter months - January 


The range of predictability depends upon the variable predicted, being 
generally shorter,' for example, for surface pressure than for 500-mb 
height. 




iKAlDil 



of 1973, 1974, and 1975 - using the so-called GISS^ GCM (Somerville 
et al., 1974). Druyan et al. (1975) have already reported evidence of 
some skill in a set of two-week averaged forecasts with this model, 
and the present study may thus be regarded as an extension of their work. 

Experiments in monthly prediction lie somewhere in a nebulous re- 
gion betoeen climate simulation and practical weather forecasting. One 
does not expect to predict successfully the synoptic details and events 
that constitute the month's weather, but one does hope that time-averaged 
properties will be at least realistically similated. Primarily, these 
experiments are designed to investigate the degree to which observed 
gross differences between the January of one year and another are deter- 
mined by the initial conditions at the beginning of each month and are 
simulated by the model, 

A description of the GISS model may be found in Somerville et al,, 
(1974) . It is a global, spherical coordinate (4 degrees of latitude 

g 

by 5 degrees of longitude), 9-layer, primitive equation, "sigma" co- 
ordinate system with upper boundary at 10 mb, integrated in 5-minute 
time steps using Arakawa's (1972) numerical method , and employing 
Arakawa's (1969) cumulus convection parameterization, parameterized solar 
and terrestrial radiative fluxes dependent on predicted water vapor and 
cloud distributions, specified surface boundary conditions^ but predicted 


The GISS model was developed at the Goddard Institute for Space Studies, 
NASA, located in New York City. 


In the version of the model used in this study the numerical method, was 
altered to include a "split grid" (Halem and Russel, 1974) in which the 
longitudinal interval increases discontinuous ly at higher latitudes . ‘At 
the same time, Arakawa's TASU (time-alternating, space-uncentered) 
scheme was eliminated. 

r 

Since. the publication of Soramerville et al, (1974), the surface albedo 
in the model has been modified to represent a greater variety of sur- 
face conditions. The climatological January albedo distribution used 
in the present study was taken from Schutz and Gates (1972) and is 
based on the original data of Posey and Clapp (1964). 
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surface land temperatures, and parameterized surface fluxes of heat, 
water vapor, and momentum, A mean state computed by the model for 
January 1973 from initial data for 20 December 1972 has already been 
published in Somerville et al. (1974) where it is compared with the 
January 5-year climatology of Oort and Rasmussen (1971). The mean 
states computed for the present study were all initialized on 1 January 
(00 GMT), and each is verified against the observed mean January state 
for the same year. 

As shown by Somerville et al. (1974), the model produces a real- 
istic, albeit imperfect, s'imulation of the mean January troposphere. 

The model stratosphere is less satisfactory, probably due to the poor 
vertical resolution at high levels. Even in the troposphere the model 
exhibits certain defects: the winter polar regions are too cold, 

meridional circulations are too weak, westerly jets are too broad, hori- 
zontal gradients are too weak, and eddy energies are too small. The 
model is still in the process of development. However, even at this 
stage , it appears to generate a sufficiently credible meteorological 
history to justify a mean monthly forecast experiment, and it is ap- 
parently representative of the current "state of the art" of numerical 
weather prediction (Druyan, 1974) . 

Two related experiments are described in this paper. In the first, 
a set of three monthly mean forecasts was computed from 12-hourly out- 
puts of the GISS model for January of 1973, 1974 and 1975, All three 
forecasts used the same fixed sea-surface temperature (SST) field, 
corresponding to the climatological SST field for January (from Washington 
and Thiel, 1970), but each was, of course, started from different initial 
atmospheric conditions. These three forecasts are discussed in sec- 
tion 2, In the second experiment, described in section 3, boundary con- 
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ditions, rather than initial conditions V7ere varied. The latter ex- 
periment, which was designed to evaluate the effect of temporal SST 
variations on the monthly mean forecasts, \ms, however, carried out 
only for January of 1974. In this study, two parallel forecasts for 
January 1974 were computed from the same initial atmospheric condi- 
tions, hut with different SST fields, one forecast using the climatolo- 
gical field and the other an observed SST field that was updated daily. 

The global data sets used for both initializ action and verifica- 
tion of the forecasts were provided by the National Meteorological 
Center (K^^G) and interpolated into the GISS grid. The January 1973 
and 1974 data were derived from a combination of the then-operational 
NMG analysis north of latitude 18N and a Flattery spectral Hough 
function analysis (Flattery,. 1971; National Weather Service, 1974) 
over the rest of the globe. For the January 1975 data, however, the 
now-operational Flattery spectral analysis was used over the entire 
globe. 

After the experiment was completed, it was learned from NMG that 

8 

an error had been committed in all the 12 GlyjT analyses for January 1975. 
Therefore, these data were subsequently discarded, and all the Janu- 
ary 1975 observed means were recomputed solely from the 24-hour ly 00 GMT 
NMG analyses. Thus, the "observed” January 1975 mean values and energetics 
were computed from 00 GMT data only, while the forecast mean values for 
all three months were computed by^ 12 -hourly averaging of both 00 and 
12_ GMT outputs. However, in order to assess the effect of 24-hourly 
vs, 12-hourly averaging, the forecast mean values for January 1975 
were computed both ways, and are discussed below. 


In the "first guess" for all 12 GMT analyses during this period, data 
for 6 October 1974 were eri'cneously accessed. This error was not 
corrected until 4 February 1975, according to NtiG staff (personal 
communication) , 
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Another, pres;imably minor, inconsistency also was committed 
in the SST update experiment for January 1974. In this calculation, 
the two parallel forecasts (one with the climatological SST and one 
with the daily updated SST) were computed using a slightly earlier 
and different version of the GISS model. Specifically, the model 
used in the SST update experiment employed an older and somewhat less 
satisfactory code for the long~wave radiation computations , as well 
as the original specification of surface albedo, with its smaller 
spatial variability (see footnote 7), Thus, the January 1974 fore- 
cast reported in section 3 below, is not quite the same as the 
January 1974 forecast described in section 2. However, in view of 
the fact that the SST update experiment was intended primarily to re 
veal the impact of temporal. SST variations and anomalies on the 
monthly forecasts, and considering other deficiencies of the model 
as a forecasting system, the inconsistency in the January 1974 com- 
putations between sections 2 and 3 is probably of trivial consequence 
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2„ EVALUATION OF THREE JANUARY FORECASTS 

In this section we examine the interannual variations of the ob- 
served January atmospheres, as they are represented by the NMC analyses, 
as well as the corresponding variations in the predicted January mean at- 
mospheres as computed by the GISS model. The forecast and observed mean 
January states are compared for each of the three years separately. As 
noted above, the three January forecasts were computed with the same model 
and identical boundary conditions. They therefore differ only because 
of the different initial conditions on 1 January of each year. Agree- 
ment between the forecast and observed atmospheres would signify (a) that 
the differences among the observed mean January states are determined by 
the initial conditions, and (b) that the model' correctly simulates this 
dependence. Disagreement could indicate that either one or both of these 
statements is false, or alternatively, that, perhaps due to defects in 
the observation and analysis system, the apparent differences among the 
"observed" .Januarys may not be real. 

The evaluation of the monthly forecasts is presented below in four 
sub-sections devoted to the energetics of the mean atmosphere, mean meri- 
dional and vertical profiles , forecast error statistics , and monthly mean 
synoptic maps. Although the forecasts were computed over the entire globe, 
only the results for the Northern Hemisphere are discussed in any detail. 
Some global data are presented, but our confidence in the global statistics 
is much lower than in the results for the Northern Hemisphere, where most 
of the observational data are concentrated. 

Energetics 

Because of the known deficiencies of the model stratosphere, the 
energetics of the predicted and observed atmospheres are computed only 
for the tropospheric region represented by the lower 8 layers of the model 
(i.e., up to about 120 mb), as 5.n- Somerville et al. (19J4) , The forecast 
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(F) and observed (0) energetics for each January, as well as the 3-year 

• 9 

averages, together with recent "climatological" estimates of January 
energetics for the Northern Hemisphere by Peixoto and Oort (1974) and 
Oort and Peixoto (1974), are presented in Table 1, (The climatological 
energetics were computed in the mixed space- time domain, while the fore- 
cast and observed energetics from this study shown in Table 1 were com- 
puted in the space domain. However, as noted by Oort and Peixoto (1974) 
and discussed by Tenenbaum ( 1976), the differences between the two methods 

of computation are considered to be small compared with errors from all 

5 -2 

other sources). The energies are given in units of 10 Joules meter 

■ -2 -2 -2 
(J-,m ) and the energy conversions in Watts meter (W m ) , Mean zonal 

available potential energy and zonal kinetic energy are represented by 

P - and K , respectively, while the corresponding eddy energies are denoted 
10 

by P and K_ . Energy conversions, P,7K,. , P„/K^ K /K , are 

EE M. MMEEEEM 

positive when the conversion is from the first form to the second form, 
and negative if in the opposite sense. Global energies are also sho^im 
in Table 1 for each January together with the 3-year averages. All 
the energy calculations for this study, both forecast and observed, are 
based on the spherical grid of 4 degrees of latitude. by 5 degrees of longi- 
tude, while the "climatological" values are based on an enlarged Hi-ic grid 
(Oort and Rasmusson, 1971). 


The "climatological" values are based on 5 years of data, 1959-1963, and 
are computed by integx'ating up to the 75 mb level. 

’^ddy energies include both transient and standing eddies. Thus , this 
part of the study is concerned not only x^ith the eddy structure of the ‘ 
mean maps, but also X7ith spatial and temporal variations during the month. 
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Table 1. Forecast (F) and Observed (0) January energetics. 

Units: energy, 10® energy conversion, Wm 

% and Fg are mean zonal and eddy available potential 
energies; Kj^ and Kg are mean zonal and eddy kinetic 
energies. Five-year "climatological" values are from 
Peixoto and Oort (1974) . 


(A) Northern Hemisphere 


1973 

1974 

1975 

11 

3"Year 

Average 

1959-1963 

Climatology 

Energy 

F 

0 

F 

0 

F 

0 

F 

0 



61.4 

55.5 

58.1 

54.5 

67.5 

55.9 

62.3 

55.3 

55.. 8 


9.0 

8.0 

8.2 

8.8 

6.8 

7o2 

8,0 

8.0 

10.5 

% 

9.9 

9.4 

9.5 

8.7 

10.7 

.7.8 

10.0 

8,6 

8.0 

% 

6 .4 

6.9 

6.5 

7.6 

4.7 

6.7 

5.8 

7.1 

9.3 


Conversion 



1 

“• 1 «. 1 ~h6 . X 1 

i 

-1,1 +0.4 

-1.2 +3.9 

-0.1 

V^E +3-1 +2.2 

+1 • 3 -f-3 « 0 i 

+2 o 3 +2 * 6 

+2.2 +2,6 

+2 • 8 

Fg/Kg +2.5 +4.0 

+2.3 +4.4 

+2,5 +0.8 

+2.4 +3.1 

+3.4 

Kg/Kj^ +0.5 +0.'3 

+0.1 +0,3 

-XO • 5 ~X0 • 3 

+0 , 4 +0 • 3 

•XO *3 


(B) Globe 


t ’ 

1 , - 
1973 

1974 

1975 

3 -Year 
Average 

F 0 

F 0 

F 

0 

F 0 

, Pj^ 45.1 43.1 

41.8 44.4 

47.3 

42.0 

44.7 43,2 

p 6.5 5.2 

E 

6.5 5.8 

4.9 

4.9 

6.0 5.3 

2,7 7.3 

7.2 7.4 

8,0 

6.5 

7.6 7,1 

Kg 5.4 4.5 

5.1 5.1 

3.8 

4.7 

5.1 4.8 


Due to a computer problem, the January 1975 forecast was run for 
only 29 days and the observed January 1975 results are for a 30 day 
period. 
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The 1975 forecast values listed in Table 1, like all the values for 
■1973 and 1974, were computed by 12-hourly averaging, using both 00 and 
12 GMT outputs, whereas, for reasons explained kbove, the 1975 observed 
values are based on 24~hourly averaging of the 00 GMT analyses only. 

The differences between the 12 -hourly and 24 hourly averaged forecast, 
mean energetics are, however, generally insignificant, as shoi-m in Table 2. 
Also sho^ra in Table 2 are the 12-hourly averaged observed mean energetics 
(including the effects of the contaminated 12 GMT analyses for which the 
wrong first guesses were used) compared with the 24-hourly averaged 
(00 GMT only) observed means. The combined effect of both the analysis 
error and the sampling intemral on the observed mean energetics is also 

* i 

seen to be relatively small. Thus, in view of the data in Table 2, 
the effect of sampling interval on the results shot-m in Table 1 may be 
ignored . 

Compared with the 5-year January climatology for the Northern Hemi- 
sphere, the average observed January data in Table 1 for the 3-year 
period 1973-1975 appear to exhibit relatively low eddy energies, which 
may or may not be the result of differences in methods of analysis and 
computation (Tenenbaum, 1976). (The use of space vs. space-time domain 
energetics accounts for only a trivial part of the difference between 

r • 

the two periods). The average observed zonal energies, on the other 
hand, are in fairly good agreement with the values from Peixoto and 
Oort (1974) . The observed conversion rates for the 3-year period are 
also in good agreement with the 5-year climatology, with the exception 
of the conversion, which is generally considered to be unreliable 

(Tenenbaum, 1976) . 


A more comprehensive diagnosis of the GI3S model's energetics, including 
a spectral analysis, may be found in Tenenbaum (1976). 
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Table 2. Comparison of 12-hour averaged (00 and 12 GMT) and 

24Thour averaged (00 GMT only) monthly mean energetics 
for January 1975. (See Table 1 and text for explanation 
of symbols) . 




F 


0 


12-hour 

24 -hour 

12-hour 

24-hour 

Energy 


Northern Hemisphere 


% 

67.5 

67.5 

55.5 

55.9 

"e 

6.8 

6.8 

7.2 

.7.2 


10.7 

10,7 

.7.8 

7.8 


h.l 

4.7 

6.3 

6.7 

Conversion 






-1.1 

-0.7 

+0.4 

+0 . 4 

V"e 

+2.3 

+2.3 . 

+2.5 

+2.6 


+2.5 

+2 ,4 

+0 . 6 

+0.8 


+0.5 

+0.4 

+0.3 

+0.3 

Energy 


Globe 



"m 

47.3 

47.4 

41.8 

42.0 

^E 

4.9 

4.9 

5ol 

’ 4.9 


8.0 

8.0 

6.5 

6.5 


3.8 

3.8 

4.5 

4.7 
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The average forecast zonal potential and zonal kinetic energies for 
the 3-year period are both higher than observed in the Northern Hemi- 
sphere. Forecast and observed zonal energies are in closer agreement, 
however, in the Southern (summer) Hemisphere, and the effect of this is 
reflected in the global values shown in Table 1. The 3-year average fore- 
cast eddy potential energy is in good agreement with the corresponding 
observed value in the Northern Hemisphere, but the forecast eddy kinetic 
energy is too low. On the other hand, the forecast eddy energies are 
higher than "observed" in the Southern Hemisphere, and the effect of 
this is also reflected in the tabulated global values. However, it appears 
quite likely that the low eddy energies apparently observed in the Southern 
Hemisphere may be simply the result of data paucity and excessive smoothing 
in the analysis. The average forecast conversion rates (except for 
are in reasonable agreement with the observed values. 

In the Northern (winter) Hemisphere, the model tends to overpredict 
the mean meridional temperature gradient each January, and hence also 
the mean zonal potential and kinetic energies. (For example, the aver- 
age vertically- integrated mean temperature difference between latitudes 
10 N and 70 N is forecast to be 30° K and observed to be 26° K for 
the 3-year period). As sho™ by Stone et al. (1975), this defect in 
the model is apparently due to an underestimate of the influx of sensible 
heat to high latitudes by large eddies, as a result of which the Arctic 
regions are forecast to be too cold in winter j The model also consis- 
tently underpredicts the eddy kinetic energy each year in the ■ 

Northern Hemisphere, a defect apparently common to all coars e-grid GCM's. 

The interannual variations of zonal available potential energy 
predicted by the model are much larger than observed in the Northern 
Hemisphere. As measured, for example, by their mean deviations, the 
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forecast and observed values of P , in the Jsforthern Hemisphere vary by 

5 -2 ^ 

3.4 and 0.5 (x 10 J m ), respectively, over the three Januarys. 

(The corresponding standard deviations are about twice as large.) For 

each of the other energy forms, the interannual variabilities of the 

forecast and observed values are more alike, with forecast and observed 

5 -2 

mean deviations (in units of 10 J m ) of 0.8 and 0.5 for P , 0.4 and 

E 

0.8 for Kj^, and 0.8 and 0.4 for K^, respectively. Thus, the model ap- 
pears to forecast greater interannual variability of P , P , and K than 
observed, but less variability of K^. The latter result, however, is 
due primarily to the relatively low value to which the observed fell 

in January 1975. From simple geostrophic (thermal wind) considerations, 
one would expect P^^ and to rise and fall together, as both depend on 
the meridional temperature gradient » However, in January 1975, the ob- 
served value of P in the Northern Hemisphere rose to its maximum for 
M 

the 3-year period, while that. of ^ fell to its minimum level. While 
this curious result could be an effect of the non-linear dependence of 
on wind speed, it is worth noting that the forecast values of P^ and 
do show the expected parallel variations of P^ and from year to 
year. The observed values of both and P^ in January 1975 are actually 
closer to the 5-year climatological values than are those of the prior 
two years, and therefore cannot be rejected as anomalous. Nevertheless, 
the possibility exists that errors in the estimates of P^ or either 
in 1975 or in the prior too years, could account for the lack of parallel 
variation of the observed values. 
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The forecast for January 1975 gives anomalously high values for the 
zonal energies and low values for the eddy energies (particularly kinetic) 
in the Northern Hemisphere. The observed eddy energies in January 1975 
are, in fact, low for the 3-year period (and much lower than the climato- 
logical values), although not as low as predicted. Thus, the model ap- 
pears to indicate correctly the sign of the negative eddy energy anomaly 
in 1975, but overpredicts its magnitude. 

In general, the model appears to forecast somewhat larger energy 
fluctuations than are observed from year to year over the Northern Hemi- 
sphere, It fails to reproduce consistently the observed pattern of inter- 
annual energy variations, indicating the correct sign of the interannual^ 
change for and over the 3-year period, but the wrong sign half the 
time for and P^. 

Meridional and Vertical Profiles 

Forecast (F) and observed (0) meridional profiles of the mean zonal 
circulation for each January are shovjn in Fig. 1 as dashed and solid 

curves, respectively, plotted as functions of latitude from pole to pole. 

-1 

The mean zonal winds, in meters sec , are averaged over pressure (height) 
as well as longitude at intervals of 4 degrees of latitude, and include 
all 9 layers of the model. For reference, the climatological 5-year mean 
zonal winds, taken from Oort and Hasmusson (1971), are also shown (as 
crosses) for the region between latitudes 10 S and 75 N. (Positive 
values denote westerlies, negative values easterlies). 

It is apparent from Fig. 1 that the model simulates the general 
configuration of the mean zonal wind profile rather realistically over 
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the 3 year period, particularly in 1973 and 1974. However,, in January 
1975 the model fails to predict what appears to be an anomalous ob- 
served mean zonal circulation in the Northern Hemisphere, which is 
characterized by a broad, weak wind maximum shifted poleward relative 
to the climatological profile, VHiile the maximum observed mean westerlies 
decreased markedly from January 1973 to 1975, the forecast indicates 
stronger maxiratim westerlies in 1975 than in 1973, with little shift in 
latitude. (On the other hand, the model does correctly simulate the 
strengthening of the westerlies in January 1975 in the higher middle 
latitudes of the Northern Hemisphere, and also indicates a somewhat 
broader maximum than in the previous two years.) The discrepancy in 
Fig. 1 between the forecast and observed mean zonal wind profiles over 
the Northern Hemisphere in January 1975 illustrates the source of the 
difference bet\-Een the forecast and- deserved mean zonal kinetic energies 
in 1975 noted in Table 1. - 

In the Southern Hemisphere, where the observed profiles in Fig, 1 
appear to indicate a biennial oscillation in the speed of the maximum 
mean westerlies, the model predicts little or no interannual variation, 
so that the forecast is out of phase with the observed variation. On . 
the other hand, the interannual shifts in the latitude of the maximum 
mean x^esterlies in the Southern Hemisphere are simulated rather well by 
the model. 

As a further illustration of the model's simulation of interannual 
variations in the large scale circulation. Fig, 2 shows the forecast 
(F) and observed (0) vertical profiles of the mean zonal wind for 
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each- January at the latitude of the observed jet stream in the Northern 
Hemisphere (latitude 30 N in January 1973 and 1974, and 34 N in 
January 1975), Again, ;the data from Oort and Rasmusson (1971), in this 
case for latitude 30 N, are plotted (as crosses) for comparison. 

The vertical wind profile, like the meridional profile of the zonal 
wind, is realistically simulated, in general, by the model over, .the 
3-year period. However, the model fails to reproduce the marked decline 
in the jet stream velocity from January 1973 to 1975, and, in fact, 
predicts a slight strengthening of the jet. Actually, the model fore- 
casts relatively little variation in the mean zonal wind profile from 
year to year compared with the observed wind analysis, and this is, of 
course, reflected in the behavior of the zonal kinetic energy, as 
noted earlier. 

The relatively weak observed circulation over the Northern Hemi- 
sphere in Januarj' 1975 noted in Figs. 1 and 2, and indicated also by 
in Table 1, appears to be due largely to a reduction of vertical 
wind shear in the upper troposphere over the sub-tropics . As shown by 
Wagner (1975), the zonal flow at 700 mb in January 1975 was actually 
"stronger than normal" over much of the Northern Hemisphere in middle 
and high latitudes, especially in the western hemisphere. Fig. 1 
does, in fact, indicate that both observed and forecast wes.terlies were 
stronger than normal in the vicinity of latitude 50 N in January 1975. 
Hovyever, the observed sub-tropical wind maximum collapsed in 1975, and 
it is this latter feature that is poorly simulated in both the meridional 
and vertical forecast wind profiles. The possibility that the 1975 "ob- 
served" wind profiles and energetics in the Northern Hemisphere may 
not be correctly represented by the newly-adopted Flattery analysis at 
NMG cannot, of course, be discounted. 
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Interaimual variations in the observed and predicted tropospheric 
mass fields are illustrated in Fig. 3, showing the mean meridional pro- 
files of geopotential height at the 505-mb level, one of the nominal 
isobaric levels corresponding to the model's sigma- coordinate system. 
Again the general configuration of the mid-tropospheric isobaric surface 
is realistically simulated. The apparent discrepancies between forecast 
and obseirved heights in the Antarctic are probably due in part to the 
model's generation of an excessively warm Antarctic in summer and in 
part to a defect in the NMC analysis over that region. The anomalously 
high elevation of the isobaric surface' observed in the Arctic in Janu- 
ary 1974 is not predicted by the model, but othert-7ise the profiles are 

. \ 

in good agreement. 

As shown above, as well as in Somerville et al. (1974), the model 
does produce a realistic mean meridional structure for the month of 
January, although it fails to simulate realistically certain interannual - 
variations of the zonally averaged global atmosphere. A more critical 
test of the model, however, would be of how well it reproduces the 
synoptic structure of the atmosphere as represented, for example, by 
the eddy energies. In F.ig. 4 are shoT-m the three January meridional 
profiles, for the Northern Hemisphere only, of the observed and pre- 
dicted eddy kinetic energies averaged over pressure (height) and longi- 

5 “2 -1 

tude. The units are 10 J m bar . As noted earlier, the eddy kinetic 
energy includes both transient and stationary eddies, so that this quan- 
tity reflects the temporal as well as spatial variability of the atmo- 
sphere x%^ithin the month, and not merely the eddy structure of the mean 
maps. 
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The interaimual variations in eddy kinetic energy noted in Table 1 
are now seen in mcnre detail in the observed energy profiles (solid 
curves) in Fig. 4. January 1974 was the most active of the three months 
in terms of (but not 
profile in the vicinity of latitudes 45-50 N, about 10 degrees north of 
the maxima in the profiles for the other two years . (The 1974 profile 
is actually more similar to the 5-year distribution of Eeixoto and 

Oort (1974) than are the curves for 1973 and- 1975.) However, the model 
fails to reproduce this peak in the energy profile, and indeed forecasts 
a higher mid- latitude maximum in 1973 than in 1974. In general, the 
model underpredicts the eddy kinetic energy in middle latitudes of the 
Norther.n_Kemisphere each January. The discrepancy betoeen observed and 
forecast eddy kinetic energies is most apparent in January 1975 when the 
model underpredicts the eddy kinetic energy every^^here in the northern 
Hemisphere outside the tropics to, yield the very loxj minimum shown 
in Table 1, In the tropics, on the other hand, the model generates too 
much eddy kinetic energy each January compared with the observed atmo- 
sphere. 

In summary, the model does not yet appear to be capable of account- 
ing for the observed interannual variations in the mean meridional 
structure and circulation of the mean January atmosphere on the basis of 
the initial state of the atmosphere at the beginning of each month. 

Error Statistics 

The monthly mean forecasts have been evaluated in terms of rms 
errors and si skill scores (Teweles and Wobus, 1954), the latter being 
a dimensionless measure of the difference between predicted and observed 
horizontal gradients. (As in the case of rms errors, lower skill scores 


K^) , and a high maximum is found in the 1974 
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signify better forecasts.) Both rms errors and SI scores were computed 

for the sea-level pressure and 500-mb height fields, but only the rms 

errors were computed for the 850-mb temperature forecasts. 

Verification statistics for each January forecast (F) are shoxvn 

in Tables 3, 4, and 5 for seven differenc regions; the whole globe, the 

Northern Hemisphere, the tropical belt (22 N to 22 S), an eastern Pacific- 

13 

United States region between latitudes 30-54 N and longitudes 7.5 W - 180^ 
North America between 30-70 N and 75-130 W, the United States between 
30-54 N and 75-130 W, and a European region between 34-86 N and 10 W - 
40 E. For the sea-level pressure and 500-mb height verifications, only 
land points were used in the North American, European, and United States 
verification regions, ’ ' 

For comparison, rms errors and SI scores were also computed for 
two "no skill" forecasts represented by climatology (M) and persistence (P). 
(The global climatology "forecast", data for mean January were provided 
by the National Center for Atmospheric Research.) The persistence "fore- 
cast" for each January was taken to be the initial state of the atmosphere 
on 1 January (00 GOT) of that year. 

In Tables 3, 4, and 5, the "best" forecast is underlined for each 
verification statistic, each forecast element, and each region. The model 
may be considered to be "successful" in a predictive sense only if it 
exceeds in skill. both climatology and persistence. Of the 87 forecast 
verifications listed in the three tables, 53, or 617o, may be rated as 
"successful" by this criterion. Of course, the 87 forecasts are not in- 
dependent. If they were, one could expect one-third to be "successful" 
just by chance on the assumption that F, M. and P are equally likely to 
excel, and one would have to concede some skill to the model. A more 

I • 

13 ~~ ' 

SI scores were not computed for this region. 
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Table 3. (A) Root-mean- square (RMS)' errors and (B) SI skill scores 

of forecast (F) mean sea-level pressure (mb) for January 
1973, 1974j and 1975, M and P denote climatology and 
persistence "forecasts", respectively. Minimum values are 
underlined. 


A, RMS Error (mb) 


Region 

1973 

1974 

1975 

■ ' 

F 

M 

P 

F 

M 

P 

F 


P 

Globe 

8.5 

10.9 

9.2 

7.6 

7.6 

9.2 

6.0 

5.9 

5.5 

Northern Hemisphere 

10.0 

8.7 

12.2 

8.6 

9.2 

11.7 

6.1 

6,6 

7.2 

Tropics 

A. 6 

8.7 

4.3 

3.7 

3.2 

1.9 

3.1 

3.0 

1.6 

E. Pacific-U.S. 

15.3 

6.2 

15.0 

8.9 

9.3 

12,4 

5.8 

6.8 

7.1 

North America 

12.8 

7.1 

12.5 

4.7 

5.5 

10.8 

7.2 

10.0 

12.4 

United' States 

10.1 

6.7 

10o6 

4.4 

3.4 

9.5 

5.6 

7.1 

8.2 

Europe 

4.2 

9.6 

15.3 

11.0 

15.5 

10.1 

7.7 

4.1 

7.3 


B. Si Score 


Region 


1973 

1 

1974 

1975 


F 

M 

P ^ 

. F 

M 

P 

F 

M 

P 

Globe 

83 

80 

73 

76 

80 

74 

74 

72 


Northern Hemisphere 

81 

81 

77 

21 

89 

81 

72 

73 

21 

Tropics 

77 

75 

68 

71 

80 

60 

71 

67 

50 

North America 

96 

95 

91 

92 

106 

90 

92 

97 

91 

United States . 

97 

98 


21 

101 

97 

95 

103 

21 

Europe 

73 

87 

69 


110 

95 

77 

55 

76 
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Table 4. (A) E.oot-mean-square (RMS) errors and (B) Si skill scores 

of forecast (F) mean 5.00-mb geopotential height (m) for 
January 1973, 1974, and 1975, M and P denote climatology 
and persistence "forecasts", respectively. Minimum values 
are underlined. 


A, RMS Error (m) 


Region 

1973 

1974 

1975 

■*. . '* * 

F 

M 

P 

F 

M 

P 

F 

M 

P 

ciobe ■ *'■ 

73 

97 

89 

2§ 

88 

93 

64 

68 

EE 

Northern Hemisphere 

72 

94 

119 


108 

116 

63 

82 

73 

Tropics 

36 

90 

35 

13 

25 

29 

IE 

34 

23 

E. Pacific-U.S. 


118 

175 

§8 

103 

114 


123 

90 

North .-America 

IZ 

132' 

146 

84 

S3 

151. 


130 

117 

United States 

% 

142 

127 

79 

80 

101 


117 

108 

Europe 


80 

169 

98 

252 

121 

113 


92 


B. Si Score 



1973 1 

1974 

1975 


F 

M 

P 

F 

M 

P 



P 

Globe 

50 

54 

55 

51 

55 

58 


49 

47 

Northern Hemisphere 


55 

62 

- ^ 

60 

64 

42 

52 

53 

...Tropics 

12 

73 

70 

62 

72 

69 

IZ 

71 

58 

ViNoxth .^America 

34 

55 

64 

38 

43 

57 

34 

50 

54 

I'-United iS t ate s 


48 

65 

35 

41 

52 

32 

40 

50 

-/'Europe - 

46 

58 

69 

59 

84 

79 

53 

36 

68 
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Table 5, Root-mean-square errors of forecast (F) mean 850-mb 

temperature (deg. C) for January 1973, 1974, and 1975, 
M and P denote climatology and persistence "forecasts" 
respectively. Minimum values are underlined. 


Region 

1973 

1974 

1975 


F M P 

F M P 1 

F M P 

Northern Hemisphere 

1 4.1 4.3 4.6 

4.7 5.1 4.7 

4.1 4.5 3.6 

E. Pacific-U,S. 

4.3 6.5 5.2 

4.4 3.3 6.3 

4.9 6.9 5.1 

United States 

4.6 7.8 5.4 

4,6 4.0 7.8 

3.5 7,5 5.2 
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interesting result, however, is revealed, by considering separately the 
the 39 sea- level pressure, the 9 850-inb temperature, and the 39 500~mb 

height forecasts. For these* three variables, the "successful" fore- 
casts are found to total, respectively, 13 (33%), 6 (677») » and 34 (877o)» 
Thus, it appears that, by this criterion, the model exhibits no skill 
over chance in forecasting monthly mean sea-level pressure, some skill 
in forecasting monthly mean 850-mb temperature, and considerable skill 
in forecasting the monthly mean 500-mb height field. 

The average rms error of the three 500-mb height forecasts over 
the Northern Hemisphere (all "successful") is found to be 72 m for the 
model, compared with 95 m for climatology and 103 m for persistence. 

The average Si skill score of the model at 500-mb over the Northern 
Hemisphere for the three January forecasts (all "successful") is 47, 
compared with 56 for climatology and 60 for persistence. The theoretical 
range of SI skill scores extends from zero, for a perfect forecast, to 
a maximum of 200. However, National Heather Service forecasters consider 
an SI score of 20 to be "virtually perfect" and a score of 70 to re- 
present a "worthless" forecast (Shuman and Hovermale, 1968). Thus, as 
a measure of "percentage skill", the Weather Service frequently employs 
the quantity (70 - Sl)/50. In terms of this quantity, the model's skill 
in predicting the monthly mean 500-mb height field over the Northern 
Hemisphere is thus 46%., compared with 28%. for climatology and 20% for 
first day persistence, (It is noteworthy that, by the criterion above, 
none of the model’s monthly mean sea-level pressure forecasts would be 
rated better than "worthless",) 
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The model's performance is particularly outstanding in the prediction 
of the 500-nib height field over the. United States, where the mean ruis error 
for the three Januarys is 55 m (compared with 113 m for climatology and 
112 m for persistence), the mean Si score is 31 (compared with 43 for 
climatology and 56 for persistence), and the mean percentage skill is 
787o (compared with 54% for climatology and 28% for persistence.) 

Prognostic and Observed Mean Maps 

To illustrate the synoptic output of the model, maps of the mean pre- 
dicted and observed sea- level pressure and 500-mb height fields are shoxm 
below for each January. Although global maps are presented, the analyses 
in the Southern Hemisphere are of rather dubious quality and the forecasts 
there are somewhat degraded by model deficiencies in the Antarctic. There- 
fore, in the interests of brevity, and in recognition of the uncertainties 
in the Southern Hemisphere, the discussion of the synoptic results in this 
section is limited to the Korthern Hemisphere only. 

Figs. 5, 6, and 7 display the global forecast and observed monthly 
mean sea- level pressure fields for each January, with isobars drawn for 
an interval of 4 mb. Figs. 8, 9, and 10 illustrate the corresponding 
500-mb height fields , xcith geopotential height contours drawn for an in- 
terval of 100 meters. All fields have been slightly smoothed with a 
single 5-point smoothing pass. 

In January 1973 the principal defects in the sea- level pressure 
prognosis (Fig. 5) for the Northern Hemisphere are found in the Aleutian 
and Icelandic regions, X'jhere the depths of the major mean c3’^clones are 
underpredicted. The model, in fact, generates an over-developed anti- 
clone in the central North Pacific, extending from the subtropics to the 
Aleutian Islands, in place of the observed North Pacific low. On the 
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other hand, the Korth Atlantic pressure pattern is simulated fairly well, 
.except for the cyclone depth, as is the Eurasian high pressure system. 

The prognostic map is generally unsatisfactory over the Northern Hemi- 
sphere, as’ “indicated’ also by the poor error 'statisti'ds 'for January' 1973 
sho™ in Table 3. 

The prognostic sea-level pressure field for January 1974 (Nig» 6) 

is again mainly deficient in it^ simulations of the deep cyclone centers 

in the North Atlantic and North Pacific, particularly the latter, as 

well as the strong, Siberian anticyclone. The predicted pressure field 

in the western North Pacific must again be rated as a failure, v?hile that 

of the North Atlantic may be considered moderately successful, although 

\ 

neither the ab'normal depth nor the location of the Icelandic lox'f are ac- 
curately forecast. On the other hand, the troughs over Mexico and south- 
east Asia and the subtropical anticyclones are represented reasonably well 
by the model, and the North Atlantic cyclone is predicted (correctly) to 
be deeper than normal. 

In January 1975 (Fig. 7), the model once again fails to predict 
correctly the location and depth of the lox^ in the North Atlantic, as 
well as the Siberian anticyclone, on the mean sea- level pressure map. 

The result is a poor forecast of the pressure gradient over Europe- In 
the North Pacific the model predicts an overdeveloped subtropical high 
displaced too far to the x-yest, and fails to- account for the mean trough 
development in the central and xjestern Pacific. 

Despite the general inadequacy of the prognostic mean -sea level 
pre.ssure maps (reflected also in Table 3), the model does siucceed in re- 
producing certain characteristics of the year-to-year variations in the 
mean January sea-level pressurd^-^eld'. The Icelandic lox^, for example 
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was observed to be relatively weak in 1973 and 1975 and very intense' in 
1974, with interannual changes in mean central pressure of about - 16 mb 
from 1973 to 1974 and +12 mb from 1974 to 1975. The forecast sea-level 
pressure maps indicate parallel interannual changes of about -20 mb 
from 1973 to 1974 and +12 mb from 1974 to 1975. Thus, while the model 
fails to predict the absolute magnitude of the pressure system in the 
North Atlantic, it appears to be capable of simulating an important 
interannual change in, the mean January circulation pattern at sea level. 

In contrast to the sea-level pressure fields, the monthly mean 
500-mb height fields are -predicted rather realistically by the model in 
all three months. Comparisons of the forecast and observed 500-mb con- 
tour patterns (Figs. 8, 9, 10) clearly illustrate what has already been 
shown in the error statistics of Table 4. Despite some obvious defects, 
such as phase shifts in the long waves, the agreement between forecast 
and observed fields is generally s'atisfactory. Furthermore, as indicated 
by Table 4, the model is not merely forecasting persistence or climatology, 
for, in terms of both rms errors and skill scores, the forecasts are 
superior to both climatology and persistence. The apparent skill of the 
model in the prediction of monthly mean 500-mb geopotential height fields 
is an encouraging result, which may have some practical value for long- 
range weather forecasting. 

Of course, the 500-mb forecasts are not without serious defects. 

In January 1973 (Fig. 8), for example, the predicted wavelength over North 
America is too long, and the ridge forecast to lie over Alaska is actually 
found much farther east on the observed map. For January 1974 (Fig. 9) 
the strength of the northern hemisphere westerlies appears to be under- 
predicted. Furthermore, the 1974 vjesterlies are displaced too far south 
in the North Atlantic and too far north in the North Pacific on the pre- 
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dieted map compared with the observed field, the amplitude of the pre- 
dicted trough over the Black Sea region is too flat, and the closed 
circulation in the western Pacific is not predicted at all. On the 
other hand, the broad, flat, zonal current over the western hemisphere 
in Januari"^ 1974, extending from the west coast of North America to the 
east coast of Europe, is simulated fairly well in the prognostic mean 
500-mb map for 1974. The January 1975 500-mb forecast (Fig, 10) for 

the Northern Hemisphere is the most successful of the three in terms 
of both the error statistics in Table 4 and the synoptic mean maps. 
Phase agreement between the predicted and observed mean contours is good 

in middle latitudes, but less satisfactory in low and high latitudes. 

\ 

Major defects are found, however, north of Japan, where the trough 
amplitude is underpredicted, and east of Greenland, where the ridge 
amplitude is overpredicted. 

In summary, the model generates a credible monthly mean prognostic 
map for the 500-mb level, and exhibits general predictive skill in ex- 
cess of climatology and persistence, despite some obvious errors. At 
sea level the model is less successful, showing no general skill in 
forecasting the monthly mean sea-level pressure field, although it did 
succeed in reproducing the interannual changes in the depth of the mean 
January Icelandic low from 1973 to 1975. 
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3. SST UPDATE EXPERIMENT: January 1974 

One of the many causes of forecast error in numerical weather pre- 
diction is incorrect specification of the sea- surface temperature (SST) 
field. '['Hiether climatological or initial observed ocean temperatures are 
used for the calculation of surface fluxes over water, the SST field will 
eventually be in error unless it can be predicted. A successful coupled 
air-sea model may help to solve this problem. On the other hand, it is 
possible that an interactive model, through positive feedback of error, 
may have even worse predictability characteristics than one in which the 
SST field is constrained by climatology. Furthermore, even if a perfect 
SST forecast could be made, the influence of SST variations over the 
forecast period may prove to be negligible compared with other causes of 
decay of predictability (Spar and Atlas, 1975). 

Some insight into the possibility of extending the useful range of 
atmospheric predictions by means of coupled models may be gained through 
experiments with a non-interactive atmospheric model in which the SST 
field, while prescribed, is altered during the forecast run to correspond 
to the observed ocean temperatures. In such an experiment, the atmo- 
spheric forecast is computed almost as it would be with a coupled model 
in which the ocean prediction provides a perfect SST forecast for the 
atmospheric calculations, but with, the feedback simulated through the use 
of observed SST’s. 

In the GISS SST update experiment, two parallel forecasts were 
computed for January 1974 from initial data for 00 GMT, 1 January, 

For the Control (C) run the climatological January SST field was used 
for the entire month, while for the Anomaly (A) run the specfied SST 
values over the global ocean areas were updated at the beginning of each 
forecast day with the appropriate "observed" values at each gridpoint. 
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As noted earlier. The C forecast computed for January 1974 and de- 
scribed in this section os not quite the same as the January 1974 fore- 
cast ('F) referred to in section 1 above, the latter having been run with 
an improved infrared radiation code and a geographically variable con- 
tinental albedo. The F forecast computed with the "modified" model 
(see section 2) exhibits some marked improvements over the C forecast, 
notably a deeper and more realistic Icelandic low in the monthly mean 
sea- level pressure field, as well as some greater deficiencies, e.g., a 
weaker and less realistic Asiatic high pressure cell. In terms of rms 
errors and SI skill scores, however, the original C forecast is actually 
somewhat better than the January 1974 F forecast made with the modified 
model over both the globe and Northern Hemisphere, although not over all 
sub-regions, as can be seen by com.paring Tables 3-5 with Tables 7-9 
below. For the purpose of assessing the impact of daily updated SST's 
on forecast quality, therefore, the original unmodified program was used 
to compute both the C and A forecasts described in this section. 

T\-?o sets of daily SST fields were used in the experiment. One set, 
obtained from the U.S. Navy Fleet Numerical Weather Central (FNWC), is 
derived from surface ship and buoy observations, supplemented by sat- 
ellite data, and is available only for the Northern Hemisphere. The 
second set, derived from window channel infrared radiances monitored by 
meteorological satellite scanning radiometers, is available for the whole 
earth and was furnished by the National Environmental Satellite Service 

. 14 

(NESS) of NOAA , The two. SST analyses are not entirely independent, 
as both use surface and satellite data, as well as climatological in- 
formation, in the data processing. Nevertheless, the fields differ 
somewhat, and, as they complement each other geographically, both were 
used to derive a single global SST field for. each day of the month. 

“ IP , , , 1 in • 

14 

National Oceanic and Atmospheric Administration. 
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The two fields were combined by a method designed to maximize the 
obseirved SST anomaly, i,e., the deviation from climatology, which in 
this case- is the mean January SST, At the same time, excessive anomalies 
were viewed as being probably erroneous. In the Southern Hemisphere, 

NESS values were used exclusively. In the Northern Hemisphere, where 
both ENWC and NESS values were available, the value corresponding to the 
greater absolute anomaly was accepted. However, in both hemispheres, if 
the NESS value indicated an aboslute anomaly in excess of 6*^C, it was 
discarded. If neither the ENWC nor NESS value was available at a grid- 
point, the January SST climatology was used. 

In view of the well-known errors in sea-surface temperature measure- 
ments by ships (Saur, 1963), as well as the errors in SST's deduced from 
clear sky infrared radiances (Smith, et al., 1974; Nark, et al., 1974), 
it is probably not unreasonable to assign an uncertainty of ±l°c to both 
sets of values. Thus, daily SST' anomalies smaller than ±1°G are 
almost certainly in the field noise and should be ignored. However, 
even larger anomalies are not necessarily reliable, particularly if they 
are of short duration and small scale, On the other hand, there are some 
persistent and larger scale features of the SST anomaly field which are 
more credible. These can be seen most clearly in the monthly mean SST 
anomaly fields for -January 1974. ■ 

Three mean January 1974 SST anomaly maps are shown in Figs. 11, 12, 
and 13, Fig. 11 represents the global SST anomaly pattern based only on 
the "NESS" satellite data. Fig. 12 shows the SST anomaly pattern in 
the Northern Hemisphere ' derived from the FmiC data, and Fig. 13 illustrates 
the SST anomaly field resulting from the combination of NESS and FHWC 
data, which were used in the present experiment. 


45 













































In the Southern Hemisphere, where only NESS data (Fig. 11), were 
used, the SST anomaly field exhibits a banded zonal structure, with cold 
anomalies at low and high latitudes and x-?arin anomalies in middle lati- 
tudes. The largest, and geographically most coherent SST anomalies 
are fo\md in the South Pacific Ocean. The NESS SST anomaly field in 
the Northern Hemisphere Fig. 11) shows a similar general pattern of 
positive anomalies in middle latitudes, with negative anomalies in high 
and low latitudes, although it is *not as well organized as in the Southern 
Hemisphere. The largest and most coherent warm anomalies in the Northern 
Hemisphere are found in the western Pacific, according to the NESS data 
(Fig. 11). 

The mean January 1974 FM'JC SST anomaly field (Fig. 12) is seen to 
be rather different from the NESS field in the Northern Hemisphere. Major 
differences between the two are found in high latitudes in both the . 
Atlantic and Pacific Oceans, and off the east coast of North America. 

Near the Aleutian Islands and adjacent to the east coast of the United 
States, the FNVJC field indicates warmer water than does the NESS field, 
while between Newfoundland and Greenland and north of Iceland the FNWC . 
data show much colder water. Differences are also found off the west 
African coast, south of Iceland, in the mid-Atlantic, in the Gulf of 
Alaska, and in the sub-tropical Pacific. On the other hand, the princi- 
pal warm SST anomalies on the western sides of the Atlantic and Pacific 
Oceans do appear on both maps . 

The composite SST anomaly field (Fig, 13), resulting from the blend- 
ing of the daily NESS and FM'JC data for January 1974, is essentially the 
same as the NESS field in the Southern Hemisphere. In the Northern 
Hemisphere, on the other hand, the warm anomalies on the western sides 
of the oceans are (as might be expected from the blending method) larger 
both in magnitude and geographical extent on the composite map (Fig. 13) 
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than on the NESS map (Fig. 11). Thus, the composite SST anomalies, 
especially in the Pacific, exhibit even more clearly the zonal pattern 
of colder than niormal sea temperatures in the equatorial region, and 
warmer than normal sea temperatures in middle latitudes of both hemi- 
spheres. (A possible consequence of this pattern of SST anomalies 
would be a weakening of both the thermally driven Hadley cell and the 
meridional transports by the tropical mean circulation (Bjerknes, 1966). 
The model did indeed generate such a response, but the effect was small 
and could not be verified against observations because of difficulties 
in estimating the observed mean meridional winds. ) 

Daily maps of SST anomalies for January 1974 (not shox-m) derived 
from NESS and FKWG data exhibit marked differences in the Northern Hemi- 
sphere. For example, large and persistent warm anomalies are found on 
the NESS maps in the western Atlantic and Pacific Oceans, off the east 
coasts of North America and Asia, -and also in the Central Pacific. How- 
ever, the east coastal anomalies on the daily FMC maps are considerably 
smaller, weaker, and less persistent than those found in the NESS data, 
while in the central Pacific the x^arm anomaly is even larger and stronger 
in the FNWG data. 

None of the anomalies persists for the full month x<rithout change; all 
parts of the anomaly field, x(fhether in NESS or FNl-JG data, exhibit marked 
fluctuations during the month in both hemispheres. In the Southern Hemi- 
sphere, the anomaly field is initially irregular, small scale, and weak, 
then grox^s into a x^ell-organized, broad-scale system toxrards the end of 
the month. In the Northern Hemisphere, on the other hand, initallly 
strong positive anomalies in the x<restern Atlantic and in the X'jestern 
and central Pacific x-7eaken during the month. Thus , the mean fields shoxra 
in Figs. 11, 12, and 13 do not represent constant features of the sea- 
surface temperature field during January 1974. The most consistent fea- 
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ture of the January 1974 SST anomaly field, one which is found almost 
every day in some form in both the FK^'?C and NESS data, is the anomalous 
warm v;ater off the chst coast of Asia, 

One must, at this time, view the SST fields, particularly the daily 
patterns, with some skepticism. Both the observational methods and the 
techniques of analysis are imperfect, and there are undoubtedly real 
fluctuations in ocean temperatures on all scales which may or may not 
be represented in the course mesh data. It should also be noted that the 
month selected for this experiment was not characterized by unusually 
large and persistent SST anomalies, such as, for example, the 1968 anomaly 
studied by Namias (1971), Thus, we should be careful not to draw too 
general or sweeping conclusions regarding atmospheric response to sea 
temperature variations from this one experiment. 

As expected, the daily forecast skill of the model degrades rapidly 
with time regardless of the SST field used. Fig. 14 illustrates the 
growth of rms errors in predicted sea- level pressure over the Northern 
Hemisphere at 24-hour intervals for a period of 3 weeks for both the 
C (dashed curve) and A (solid curve) forecasts. For the first 5 days 
the rms ei'rors are virtually identical for the two forecasts. During the 
next 9 days the rras error of the A forecast is slightly smaller than that 
of the C forecast, but then it rises above the error of the G forecast 
during the third week. However, neither of the sea- level pressure fore- 
casts retains any skill over climatology beyond 5 days, and the differences 
between the daily C and A forecast errors are clearly of no practical 
significance. Similar results (not sho\m) are found also at the 500-mb 
level over the Northern Hemisphere, although at that level both fore- 
casts retain some skill over climatology up to about 10 days and the 
two sets of rms errors are almost identical for the three week period. 
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These results are consistent with those o£ Spar and Atlas (1975) who 
found that the use of observed SST values did not yield any detectable 
improvement in the quality of daily large-scale prognostic maps over a 
two -week period. 

The monthly mean sea-level pressure fields for January 1974 fore- 
cast with either the C or A model suffer from certain obvious defects 
as seen in Fig, 15 for the Northern Hemisphere. In the North Atlantic, 
the predicted Azores ~ Bermuda high and Icelandic low are too weak com- 
pared with the observed systems. Hence, the predicted pressure gradients 
in the North Atlantic are also much weaker than observed. In the North 
Pacific, the forecast Aleutian low is not only weaker than observed, but 
is displaced too far to the .east, so that again the predicted pressure 
gradient in the western Pacific is much weaker than observed, and the 
pressure field is, in fact, quite unrealistic. 

Comparing the A and C forecasts, one finds only a relatively small, 
and indeed negative, influence of the. SST anomalies on the predicted 
mean monthly sea-level pressure field. The Icelandic low and Azores - 
Bermuda high are both slightly weaker in the A than in the C forecast, 
and, hence, the predicted sea-level pressure gradients in the North At- 
lantic are even more in error in A than in C. The deep cyclone in the 
western Pacific is not well predicted in the A computation, being dis- 
placed into the eastern Pacific, resulting in an even less satisfactory 
sea- level pressure field than that forecast by C. However, the dif- 
ference (not shown) between the C and A forecast mean sea- level pressure 
fields are so small "that they are undoubtedly well below the noise level 
of the model, and are of no practical significance. (For a study of the 
noise level of a general circulation model, specifically the KCA.R model,' 
see Chervin and Schneider (1975 a,b).) Although corresponding noise level 
maps based on multiple random perturbations of the initial state are not 
yet available for the GISS model, it is reasonable to expect that they 
would resemble those computed for the NCAR model.) 
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Another qualitative test of the impact of a variable SST field on 
forecast skill is the degree to which the deviation of the atmosphere 
from its climatological normal state is predicted. In particular, it is 
of interest to know how the major "centers of action" in the sea-level 
pressure field in a given month depart from normal, and whether these 
departures are better predicted .by a forecast computed with a variable 
SST field than by one based on climatological SST's. In Table 6 are 
listed the latitudes, longitudes and central pressures of the five major 
sea- level pressure systems in the Northern Hemisphere, including the 
normal January positions and pressures (estimated from Crutcher and 
Meserve, 1970), the observed January 1974 values, the values predicted 
by the C and A runs, and, for comparison, the values for January 1973. 
Although they are of no particular statistical significance, the data 
sho^m in Table 6 do indicate the relative impact of the SST anomalies 
on the monthly mean forecast sea- level pressures.. 

The Icelandic low was much deeper than normal on 1974, but close to 
its normal position. Both its location and the sign of its sea-level 
pressure deviation from normal were, in fact, correctly predicted by the 
C forecast, although the depth of the Icelandic low was not. The A 
forecast, on the other hand, did not improve on the C forecast either in 
location or central pressure. The subtropical Atlantic high pressure 
belt was also close to its normal pressure and latitude in January 1974. 
^flhile both the G and A computations indicated approximately the correct 
latitude for the system, neither forecast the pressure correctl^^, and, 
of the two, the A forecast was the poorer with, regard to the deviation 
from normal. The continental anticyclone over Siberia in January 1974 
was in a nearly normal state of development, hut split into two centers. 
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Table 6, Locations and central pressures of the major 
centers of action in the Northern Hemisphere 


System 

Normal 

January 

Observed 
January 1974 

C-Forecast 

A-Forecast 

January 1973 


Lat. 

6 ON 


59N 

55N 

60N 

Icelandic 

Long. 

30W 


3 OH 

SOW 

4 OH 

Low 

Pressure 

(mb) 

996 


992 

1 

994 

986 

Azores 

Lat. 

25 - 35N 

SON 

28N 

28N 

25 - 35N 

Bermuda 

Long. 

- 

- 

- 

- 

- 

High 

Pressure 

1024 

1025 

1018 

1016 

1024 

Asiatic 

Lat. 

SON 

SON + 65N 

45N 

45W 

45N 

High 

Long. 

95E 

95E 

120E 

120E 

IlOE 


Pressure 

1034 

1032 

1024 

1024 

1030 


Lat. 

SON 

45N 

45N 


55N 

Pacific 

Long . . 

165W+ 170E 

170E 

180 


145W + 17 OE 

Low 

1 

Pressure 

i 

998 

994 

1010 ■ 


1002 

East 

1 

i 

Lat. 

3 ON 

SON 

. 30N 

3 ON 

28N 

Pacific 

Long. 

140W 

13 OW 

13 OW 

13 OW 

130W 

High 

Pressure 

1022 

1018 

1018 

1018 

i 

i 

1022 
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Both forecasts placed the high center east of its observed (and normal) 
position, and both were equally in error in the central pressure. In the 
Pacific Ocean, the western lobe of the Aleutian low was dominant in 
January 1974, and slightly deeper than normal. Both foreaasts failed to 
reflect this development, indicating a weaker than normal Aleutian low, 
displaced too far to the east, with the larger position error in the A 
forecast. The east Pacific high pressure cell on the other hand was 
equally well-predicted by both the C and A runs. Thus, in a subjective, 
qualitative sense, the use of daily updated SST's did not improve the 
prognostic monthly mean sea-level pressure field for January 1974. 

The forecast mean 500-mb height fields for January 1974 with and 
without SST updating are illustrated in Pig. 16. As in the case of the 
sea- level pressure fields, the differences between the C (top) and 
A (middle) forecast fields are relatively small compared with the dif- 
ferences between either of the forecasts and the "observed” (bottom) 
field. Both the C and A forecasts are about equally satisfactory in the 
western hemisphere and about equally unsatisfactory in the western North 
Pacific. (Compare also the forecast and observed 500-mb fields in Fig. 9.) 
Thus, again in a subjective, qualitative sense, the use of daily updated 
SST’s did not improve the prognostic monthly mean 500-mb height field 
for January 1974. 

A more quantitative demonstration of the impact of the updated 
SST's on forecast quality is presented in Tables 7-9, showing the rms 
errors and Si skill scores, respectively, for sea- level pressure, 500-irib 
height, and S50-mb temperature (rms error only) for both the control (C) 
and anomaly (A) forecasts. Also shown again for comparison are the 
climatology (M) and persistence (P) forecast scores, as in Tables 3-5. 

The geographical regions are the same as in Tables 3-5, and minimum values 
are again underlined. 
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Table 7 


. (A) Root -me an- square (RMS) errors and (B) SI skill scores of 

control (C) and anomaly (A) forecast mean sea-level pressure 
(mb) for January 1974. M and P denote climatology and 
persistence "forecasts", respectively. Minimum values are 
underlined. 


A. RMS Error (mb) 


Region 


Forecast 



C 

A 

M 

P 

Globe 

6.9 

6.9 

7.6 

9.2 

Northern Hemisphere 

7.3 

7.4 

9.2 

11.7 

Tropics 

3.6 

3.6 

3.2 

1.9 

East Pacific-U.S. 

6.6 

6,6 

9.3 

12,4 

North America 

7.5 

5.7 

5.5 

10,8 

U.S. 

5.5 

4.6' 

3.4 

9.5 

Europe 

8.3 

6.8 

15.5 

10.1 


B. SI Score 



C 

A 

M 

P 

Globe 

69 

71 

80 

74 

Northern Hemisphere 

iZ 

70 

89 

81 

Tropics 

§1 

70 

80 

60 

North America 

§9 

98 

106 

90 

U.S. 

93 

99 

101 

97 

Europe 

70 

69 

110 

95 
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Table 8. 


(A) Root-mean- square (RMS) errors and (B) SI, skill scores of 
control (C) and anomaly (A) forecast mean 500-mb geopotential 
height (m) for January 1974. M and P denote climatology and 
persistence "forecasts", respectively. Minimum values are 
underlined. 


A. RMS Error (m) 


Region 


Forecast 




C 

A 

M 

P 

Globe 


74 

88 

93 

Northern Hemisphere 

Ik 

83 

108 

116 

Tropics 

39 

40 

25 

29 

E. Pacific-U.S. 

83 

81 

103 

114 

North America 

90 

85 

83 

151 ■ 

U.S. 

100 

. 91 

80 

101 

Europe 

95 

92 

252 

121 

B. SI Score 






C 

A 

M 

P 

Globe 


49 

55 

58 

Northern Hemisphere 


49 

60 

64 

Tropics 

65 

69' 

72 

69 

North America 

18 

35 

43 

57 

U.S. 

13 

30 

41 

52 

Europe 

iz 

59 

84 

79 
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Table 9. Root-mean- square errors of control (C) and anomaly (A) 
forecast mean 850-mb temperature for January 1974. 

M and P denote climatology and persistence ’’forecasts", 
respectively. Minimum values are underlined. 


Region Forecast 



C 

A 

M 

P 


Northern Hemisphere 

4.6 

4.7 

5,1 

4,7 


E. Pacific-U.S, 

5.6 

5.1 

3,3 

6,3 


U.S. 

6:0 

5.1 

4.0 

7.8 
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With regard to the impact of SST updating on rms error. Tables 7-9 
are somev?hat ambiguous. Over the large geographical regions (globe, 
tropical belt, and Northern Hemisphere) the C forecasts show no larger 
rms errors than the A forecasts in all three prediction variables, in- 
. -dicating no beneficial impact of the SST updating. However, over smaller 
regions ("United States", "North America" , "Europe"), the rms errors are 
■"soial-i-er for the A forecasts, suggesting some possible regional beneficial 
.-'ir^^luence of SST updating. The S scores, on the other hand, with only 

- i ' 1 

.. .a^tsivof exception (sea- level pressure over Europe) indicate that the C 
■: .xf-oreca^^ts are superior to the A forecasts, and show no clearly beneficial 
effect of SST updating on the model predictions. 

The one limited prediction experiment described above thus indicates 
that updating sea-surface temperatures did not result in any clear-cut 
. '.iinprovement in forecast quality over a period of one month, either in 

the daily or monthly mean fields. . Indeed, the impact of the SST anomalies 
on the prognoses was very slight. It must of course, be recognized that 
the SST anomalies in January 19*74 were relatively modest in scale, magni- 
tude, and persistence, and are not representative of the large, persistent, 
and broad-scale anomalies which are occasionally found over the oceans. 
However, in view of the inherent decay of predictability of atmospheric 
models, it is doubtful that a meaningful test of the impact of even very 
■large .SST anomalies, or of any other influence, can be carried out for a 
±orecast period in excess of a few days until the predictive skill of the 
.models is substantially increased. The difficulty in demonstrating 
’Istratistical -significance for the impact of SST anomalies against the back- 
ground of:Tiieteorological "noise" with present-day general circulation 
models -.ii? .well illustrated in the experiments of Chervin et al. (1976).. 
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dieted map compared with the observed field, the amplitude of the pre- 
dicted trough over the Black Sea region is too flat, and the closed 
circulation in the western Pacific is not predicted at all„ On the other 
hand, the broad, flat, zonal current over the western hemisphere in 
January 1974, extending from the west coast of Korth ^.erica to the east 
coast of Europe, is simulated fairly well in the prognostic mean 500-mb 
map for 1974. The January 1975 500-mb forecast (Pig. 10) for the fJorth- 

ern Hemisphere is the most successful of the three in terms of both the 
error statistics in Table 4 and the synoptic mean maps. Phase agree- 
ment between the predicted and observed mean contours is good in middle 
latitudes, but less satisfactory in low and high latitudes. Major defects 
are found, however, north of Japan, where the trough amplitude is under- 
predicted, and east of Greenland, where the ridge amplitude is over- 
predicted. 

In summary, the model generates a credible monthly mean prognostic 
map for the 500-irib level, and exhibits general predictive skill in ex- 
cess of cliinatologyo and persistence, despite some obvious errors. At 
sea level the model is less successful, showing no general skill in fore- 
casting the monthly mean sea- level pressure field, although it did succeed 
in reproducing the interannual changes in the depth of the mean January 
Icelandic low from 1973 to 1975, 
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4. SUMMARY AND CONCLUSIONS 

The three January forecasts with the GISS model, using a fixed 
climatological SST field, have shown that, while the model, is capable 
of simulating realistically the general structure and circulation of the 
mean troposphere,, it fails to account satisfactorily for the observed 
interannual variations in either the monthly mean energetics or zonally- 
averaged circulation of the atmosphere. Thus, it must be concluded 
that the model is not yet capable of explaining the year-to-year fluctua- 
tions in the monthly mean state of the atmosphere on the basis of the 
initial conditions at the beginning of the month. Nor indeed has it 
been shovm that the observed monthly mean state is in fact determined 
primarily by the initial conditions as defined by the large-scale analysis. 

As a long-range forecasting system, the GISS model exhibits no 
general skill in predicting the monthly mean sea- level pressure field. 
However, the model did succeed in .reproducing an observed interannual 
change in the depth of the mean January Icelandic low, which was weak in 
1973 and 1975 and abnormally strong in 1974, 

The model does appear to show consistent skill in forecasting the 
monthly mean 500-mb height field over the Northern Hemisphere, The 
three January monthly mean 500-mb forecasts are superior to climatology 
and persistence in terms of both rms errors and Si skill scores, especially 
over North America, and the monthly mean prognostic maps compare favorably 
with the observed fields. 

The use of daily updated sea-surface temperatures in the prediction 
for January 1974 produced no detectable beneficial effect on the fore- 
casts, and, indeed, the total impact of observed versus climatological 
SST's on the evolution of the predicted large-scale atmospheric fields 
over a period of one month was very slight. 
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In the case of the sea- level pressure field, the decay of predictabil- 
ity of the model takes place so rapidly that, by the time the diabatic 
influence of the SST anomalies is felt by the atmosphere, all predictability 
has been lost, and no beneficial impact on monthly mean fields is demon- 
strable. At upper levels, specifically at 500-mb, where predictability 
decays less rapidly and the prognostic monthly mean fields appear to ex- 
hibit some skill over both climatology and perisitence, the physical im- 
pact of SST anomalies in the model is apparently so slight that again it 
is not possible to demonstrate any beneficial impact of the use of updated 
sea temperatures on the predicted monthly mean fields. 
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List of Figures 


1. Meridional profiles of mean zonal wind (m sec ) averaged over 
pressure (height) and longitude for January 1973, 1974, and 1975. 

Solid curves represent observed (0) winds from the NMC analysis, 
dashed curves indicate forecast (F) winds from, the GISS model, 
and crosses denote the 5-year (1959-1963) mean winds from Oort 
and Rasitiusson (1971). (Positive values denote westerlies, nega- 
tive values easterlies) . 

2. Vertical profiles of zonally-averaged mean zonal winds for January 1973, 
1974, and 1975 (m sec" ) at the latitude of the observed Northern 
Hemisphere jet stream (30 N in 1973 and 1974, and 34 N in 1975) . 

Solid curves: observed (0). Dashed curves: forecast (F) . Crosses 

denote 5-year mean winds from Oort and Rasmusson (1971). 

3. Meridional profiles of zonal mean geopotential height (geopotential 

decameters) at the 505 mb level for January 1973, 1974, and 1975. 
Solid curves: observed (0). Dashed curves: forecast (F) . 

5 -2 

4. Meridional profiles of eddy kinetic energy (units: 10 J m ) in the 

Northern Hemisphere averaged .over pressure (height) and longitude for 
January 1973, 1974, and 1975. Solid curves: observed (0). Dashed 

curves : forecast (F) . 

5. January 1973 mean sea- level pressure. Top (F) : Forecast. Bottom (0): 

Observed. Isobars are dra™ for an interval of 4 mb. 

6. January 1974 mean sea- level pressure. Top (F) : Forecast. Bottom (0): 

Observed. 

7. January 1975 mean sea- level pressure. Top (F) : Forecast. Bottom (0): 

Observed, 

8. January 1973 mean 500-mb height. Top (F) : Forecast. Bottom (0): 

Observed, Contours are draw for an interval of 100 geopotential 

.. meters. 

9. January 1974 mean 500-mb height. Top (F) : Forecast. Bottom (0): 

Observed. 

10. January 1975 mean 500-mb height. Top (F) : Forecast, Bottom (0): 

Observed. 
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11. January 1974 sea-surface teruperature (SST) anomaly (degrees C) based 
on satellite data provided by National Environmental Satellite Ser- 
vice (NESS), NOAAo (Anomaly is the observed deviation relative to 
the January mean SST field from Washington and Thiel, 1970.) 

12. January 1974 sea-surface temperature (SST) anomaly (degrees C) in the 
Northern Hemisphere based on data provided by U.S. Navy Fleet Numeri- 
cal Weather Central, Monterey, California. 

13. January 1974 sea-surface temperature (SST) anomaly (degrees C) derived 
from merger of data in Figs. 11 and 12„ (See text for details). 

14. Growth of root-mean-square (rms) errors of predicted sea-level pressure 

(mb) over the Northern Hemisphere during first 3 weeks of January 1974. 
Dashed curve: G forecast. Solid curve: A forecast. 

15.. January 1974 mean sea-level pressure, Top (C) : predicted with clima- 
tological SST’s. Middle (A): predicted with daily updated SST'so 
Bottom (0): observed, 

16, January 1974 mean 500-mb height. Top (G) : predicted with climatolog- • 
ical SST’s. Middle (A): predicted with daily updated SST's. Bottom 
(0) : observed. 
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